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Abstract. New results on hadron physics from the Belle experiment are presented. 
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O : INTRODUCTION 

These results are obtained using various data samples from 80 fb 1 to 150 fb 1 
taken with the Belle detector [1]. We identify B candidates by two kinematic vari- 
ables: the energy difference, AE = (J^,-J?i) — Ej,, and the beam constrained mass, M^ c = 

iPi . \l^b~ (Li Pi) 2 ' where E b = y/s/2 is the beam energy and p { and E t are the momenta 
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and energies of the decay products of the B meson in the CM frame. The inclusion of 
charge conjugate modes is implicit throughout this report. 
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£ : OBSERVATION OF 0+ AND 1+ BROAD CU STATES 

A study of charmed meson production in B decays provides an opportunity to test 
predictions of Heavy Quark Effective Theory (HQET) and QCD sum rules. B decays 
to D^Tt final states are its dominant hadronic decay modes and are measured quite 
well [2]. The large data sample accumulated in the Belle experiment allows to study 
production of D meson excited states. D**s are P-wave excitations of quark-antiquark 
systems that contain one charmed and one light (u,d) quark. 

B — > D**n decays have been studied using the D + n~n~ and D* + n~n~ final 
states [3]. 

Figure 1 shows the AE distributions for the B — > D + n~n and B — > D* + n~n~ 
candidates. The following branching fractions are measured: 33{B~ — > D + tz~k~) = 
(1.02±0.04±0.15) x 10~ 3 and m{B~ ^D*+%-%-) = (1.25 ±0.08 ±0.22) x 10~ 3 , 
without any assumption about the intermediate final states. 

To study the dynamics of B — > D^tztz decays, analyses of the Dalitz plots are per- 
formed. The description of the Dalitz plot D + n n includes amplitudes of the known 
D^n^ mode, possible contributions of the processes with virtual D*°n~ and B*°n~ pro- 
duction and an intermediate D + K broad resonance structure with free mass and width. 
Figure 2(a) shows the D + K invariant mass distribution together with the resulting fit. 
A clear signal of the broad resonance with J p = + is observed which can be identify as 
the scalar Dq° state. The results of the mass, width and branching fraction products are 
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FIGURE 1. AE distributions for the B -> D+%-%- (left) and B -> D*+%-%- (right) candidates. 

TABLE 1. Branching fractions and resonance parameters for the D^ 4 %~ %~ final states. 



Mode S$(B- ^D x [D^ + n-]n-), M(D X ), F(D X ), 

10~ 4 MeV MeV 

B -^Df[D+K-]n- 3.4 ±0.3 ±0.6 ±0.4 2462±2.1 ±0.5±3.3 45.6 ± 4.4 ± 6.5 ± 1.6 

B — > Dq {D + 7t ]n~ 6.1±0.6±0.9±1.6 2308± 17± 15±28 276±21±18±60 

B D^+n-}^- 6.8±0.7±1.3±0.3 2421 ± 1.5 ±0.4 ±0.8 23.7±2.7±0.2±4.0 

B -^Df[D*+n-]n- 1.8 ±0.3 ±0.3 ±0.2 [5] [5] 

B Df[D*+ 7t-]7t- 5.0 ±0.4 ±1.0 ±0.4 2427±26±20± 15 384+^ 7 ±24±70 



presented in Table 1 . 

For the Z)* + 7i; 7T final state the fit of the density distribution is performed in four 
dimensional phase space to take into account the angles of the pion from D* decay. The 
fit function includes both known D^ , D® intermediate state contributions and a broad 
D* + K resonance with free parameters. Figure 2(b) shows the D* + n invariant mass 
distribution as well as the resulting fit. Together with the narrow resonances a clear 




FIGURE 2. Minimal D + n~ (a) and D* + n~ (b) invariant mass distributions. The points with error bars 
correspond to the B meson signal events, the hatched histogram shows the sidebands. The open histogram 
is the result of a fit while the dashed one shows the fit function without narrow resonance contribution. 



signal of the broad state is observed. The angular distribution of D*7C from this state is 
consistent with J p = 1 + , j q = 1/2. This state can be identified as a P-wave excitation of 
cu - D®. The results of the mass, width and branching fraction products are presented 
in Table 1 . 

Together with observations of the broad resonances the branching ratios of B decay 
to the modes with known D**: D®n~ and D^7t~ have been measured. Using these 
measurements the ratio of D^ branching fractions h = SS(Df — > D + K~) / £3(D'f — >■ 
D* + 7t~) = 1.9 ±0.5, consistent with the world average h = 2.3 ±0.6 [2], is obtained. 
The measured ratio R = 38(BT -> Dfir)/38(B- -> D\k~) = 0.77 ±0.15 is lower 
than the CLEO measurement 1.8 ±0.8 [4] (although the results are consistent within 
errors) but is still a factor of two larger than the factorization prediction [6]. From our 
measurement it is impossible to determine whether the non-factorized part for tensor 
and axial mesons is large, or whether higher order corrections to the leading factorized 
terms should be taken into account. 

Our measurements show that the narrow resonances compose (36 ± 6)% of the Dnn 
decays and (63 ±6)% of the D*nn decays. This result is inconsistent with the QCD sum 
rule prediction and may indicate a large contribution from a color suppressed amplitude. 



OBSERVATION OF NEW STATES D+ (2317) AND D+ (2457) 

The narrow D s n° resonance at 2317 MeV, recently observed by the BaBar collabora- 
tion [7], is naturally interpreted as a P-wave excitation of the cs system. The observa- 
tion of a nearby and narrow D*7t° resonance by the CLEO collaboration [8] supports 
this view, since the mass difference of the two observed states is consistent with the ex- 
pected hyperfine splitting for a P-wave doublet with total light-quark angular momentum 
j = 1/2 [9, 10]. The observed masses are, however, considerably lower than potential 
model predictions [1 1], and similar to those of the cu j =1/2 doublet states recently re- 
ported by Belle [3]. Measurements of the D sJ quantum numbers and branching fractions 
(particularly those for radiative decays), will play an important role in determining the 
nature of these states. 

We confirmed both resonances and measured masses for + and 1+ states to be 
(2317.2 ±0.5 ±0.9) MeV and (2456.5 ± 1.3 ± 1.3) MeV respectively [12]. We also 
report the first observation of the radiative decay D sJ (2A51) — > D s y. Figure 3 shows the 

mass difference between the D^n and candidates. The ratio '^^j^asi^d*"^) * s 
found to be 0.55 ±0.13 ± 0.08. 

We also search for D s j production in B — > DD s j decays [13]. We reconstruct D°(D~) 
mesons in the K + n~, K + 7t~n~n + and K + n~n (K + n~n~) decay channels. D+ 
mesons are reconstructed in the K*°K + and K®K + decay channels. D s j candi- 

dates are reconstructed from mesons and a n°, y, or n + n~ pair. The mass differ- 
ence M{D s j) —Mip^) is used to select D s j candidates. We use central mass values of 
2317 MeV and 2460 MeV for £>. s y(2317) and D sJ (2451) respectively and define signal 
regions within 12 MeV for the corresponding mass difference. We observe a clean signal 
for B -> DD sJ (23n)[D s 7t°] and B -> DD s j (2457) [D*Jt ]. We also observe for the first 
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FIGURE 3. M(D s 7t°) -M Ds (a), M(D*n°) - M D * (b) and M(D s y) -M Ds (c) mass-difference distribu- 
tions. The signal is described using a double Gaussian and a third-order polynomial for the background. 
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The histogram shows no structure for the Dl sidebands. 




2.3 2.4 2.5 

M(D) (GeV/c 2 ) 




-0.5 0.5 

cos(e Dsy ) 



FIGURE 4. Left: M(D sJ ) distribution for the B -> candidates: (a) D^(2317) -> Z). S 7T°, (b) 

D s /(2457) -» D*w° and (c) £>,j(2457) -> D s y. Right: the D s y(2457) -> D s yhelicity distribution. Points 
with errors represent the experimental data and curves are the results of the fits. 



time the D. v y(2457) — > D s y decay. Figure 4(left) shows the invariant mass distributions 
for these decays. The measured branching fractions are presented in Table 2. We obtain 

the ratio f^^^ff, = 0.38 ±0.11 ± 0.04, which is consistent with that from the 

continuum study. 

We also study the helicity distribution for the D s .y(2457) — > D s y decay. The helicity 
angle 6d s j is defined as the angle between the D s j(2451) momentum in the B meson 
rest frame and the D s momentum in the D s j(2457) rest frame. The 9o s y distribution in 
the data (Fig. 4(right)) is consistent with MC expectations for the 7=1 hypothesis for 
the D.vy(2457) (^ 2 /n.d.f= 5/6), and contradicts the J = 2 hypothesis (^ 2 /n.d.f .= 44/6). 
The 7 = hypothesis is already ruled out by the conservation of angular momentum and 
parity in D sJ (2457) -> D s y. 



TABLE 2. B — > Z)D s j branching fractions. 



Decay channel 



38, 10 
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Signif. 



B^DD sJ (23 17) [D,w°] 
B^DD. s/ (2317) [D*r] 
B^DD sJ (2A51) [D*nP] 
B^DD sJ {2A51) [D s y] 
B-tDDsj (2457) [D*y] 
B ^ DD sJ (2451) [D s n + %-\ 
B^DD S] (7A51) [D s 7t°] 



8.5+^ ±2.6 
2.5l?;}(<7.5) 
17.8+^| ±5.3 
6.7+};i±2.0 
2.7+}'^(<7.3) 

< 1.6 

< 1.8 
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FIGURE 5. Distributions of M(n + n l + l )—M(l + l ) for selected events in the AE-M\, C signal region 
for (a) Belle data and (b) generic BB MC events . 



OBSERVATION OF A NEW NARROW CHARMONIUM STATE IN 

B ± -> K+Tt+Ti-J/yf DECAY 

A major experimental issue for the cc charmonium particle system is the existence of as 
yet unestablished charmonium states that are expected to be below threshold for decays 
to open charm and, thus, narrow. These include the n = 1 singlet P state, the /z c (lP), 
and possibly the n = 1 singlet and triplet spin-2 D states, i.e. the J PC = 2~ + 1 D C 2 and 
J PC = 2— l 3 D c2 , all of which are narrow if their masses are below the DD* threshold. 
The observation of these states and the determination of their masses would provide 
useful information about the spin dependence of the charmonium potential. 

We report on an experimental study of the k + k~J /\jf and JXcQ mass spectra from 
exclusive B + — ► K + 7t + x~J/y/ and K + yx c o decays [15] using a 152M BB event sample. 
For the B — > Kk + k J '/ ty/ study we use events that have a pair of well identified oppo- 
sitely charged electrons or muons with an invariant mass in the range 3.077 < M/+;- < 
3.1 17 GeV, a loosely identified charged kaon and a pair of oppositely charged pions. 

Figure 5(a) shows the distribution of AM = M(tz + tz~1 + 1~) —M(l + l~) for events in 
the AE-M hc signal region. Here a large peak corresponding to \j/(2S) — > Tt + n~J /*¥ 
is evident at 0.589 GeV. In addition, there is a significant spike in the distribution at 
0.775 GeV. Figure 5(b) shows the same distribution for a large sample of generic BB 



TABLE 3. Results of the fits to the y/(2S) and M = 3872 MeV 
regions. The errors are statistical only. 





Quantity 


y(2S) region 


M = 3872 MeV region 


Signal events 


489 ±23 
3685.5 ±0.2 MeV 
3.3 ±0.2 MeV 


35.7±6.8 
3871.5 ±0.6 MeV 
2.5 ±0.5 MeV 




FIGURE 6. Signal-band projections of (a) M bc , (b) M K+7[ -j/ ¥ and (c) A£ for the X(3872) -> K+K-J/y 
signal region with the results of the unbinned fit superimposed. 



Monte Carlo (MC) events. Except for the prominent ty/(2S) peak, the distribution is 
smooth and featureless. 

We make separate fits to the data in the y(2S) (3580 MeV < M n+n - J/v < 3780 MeV) 
and the M = 3872 MeV (3770 MeV < M k+jz -ji w < 3970 MeV) regions using a simul- 
taneous unbinned maximum likelihood fit to the M^ c , AE, and M n + n - } jy distributions. 
The results of the fits are presented in Table 3. Figures 6(a), (b) and (c) show the Mb c , 
Mx+x-j/y, and AE signal-band projections for the M = 3872 MeV signal region, re- 
spectively. The superimposed curves indicate the results of the fit. There are clear peaks 
with consistent yields in all three quantities. The signal yield of 35.7 ±6.8 events has a 
statistical significance of 10.3a. In the following we refer to this as the X(3872). 

We determine the mass of the signal peak relative to the well measured \y(2S) mass: 
M x = Mf eas -M^L + = 3872.0 ± 0.6 ± 0.5 MeV. Since we use the precisely 

known value of the y(2S) mass [2] as a reference, the systematic error is small. The 
measured width of the X (3872) peak is o = 2.5 ±0.5 MeV, which is consistent with the 
MC-determined resolution and the value obtained from the fit to the yf(2S) signal. From 
this we infer a 90% confidence level (CL) upper limit of T < 2.3 MeV. 

We determine a ratio of product branching fractions for B + — > £ + X(3872) [n + n~J/ \jf] 
andfl+ ^K + \j/(2S)\K + 7t-J/\j/] to be 0.063 ±0.012 ±0.007. 

The decay of the D c 2 charmonium state to YXcO lS an allowed El transition with a 
partial width that is expected to be substantially larger than that for the n + iz~ J j\\f final 
state; e.g. the authors of Ref. [16] predict r( 3 D c2 — > YXco) > 5 x r( 3 D c2 — ► 7t + 7t~ J/yr). 
Thus, a measurement of the width for this decay channel can provide important infor- 
mation about the nature of the observed state. We searched for an X(3872) signal in the 
YXco decay channel, concentrating on the XcO ~^ jJ/ ¥ fi na l state. 

We select events with the same — > l + l~ and charged kaon requirements plus 
two photons, each with energy more than 40 MeV The signal-band projections of M^ c 




and My Xc0 for the yf(2S) region are shown in Figs. 7 (a) and (b), respectively, together 
with curves that represent the results of the fit. The fitted signal yield is 34. 1 ± 6.9 ±4.1 
events. The number of observed events is consistent with the expected yield of 26 ±4 
events based on the known B — > K\ir(2S) and ty(2S) — > YXcO branching fractions [2] and 
the MC-determined acceptance. 

The results of the application of the same procedure to the M = 3872 MeV region 
are shown in Figs. 7(c) and (d). Here, no signal is evident; the fitted signal yield is 
3.7±3.7±2.2. From these results, we determine a 90% CL upper limit on the ratio 
of partial widths of njp^p^rg^ < 0.89. This limit on the fXcO decay width 

contradicts expectations for the 3 D C 2 charmonium state. 

The mass of the observed state is higher than potential model expectations for the 
center-of-gravity (cog) of the l 3 D cJ states: M cog (lD) = 3810 MeV [19, 18]. 

In summary, we have observed a strong signal for a state that decays to n + n~J/\if 
with M = 3872.0 ±0.6 ±0.5 MeV and T < 2.3 MeV (90% CL). This mass value and 
the absence of a strong signal in the yXcO decay channel are in some disagreement with 
potential model expectations for the 3 D C 2 charmonium state. The mass is within errors 
of the D°D*° mass threshold (3871.3 ±1.0 MeV [2]), which is suggestive of a loosely 
bound DD* multiquark "molecular state," as proposed by some authors [17]. 



MEASUREMENT OF THE E + E~ -> £>(*)+£>(*)- CROSS-SECTIONS 

The processes — > D^ + D^~ have not previously been observed at energies 
y/s ^> 2Mr>. A calculation in the HQET approach based on the heavy-quark spin symme- 
try [20], predicts cross-sections of about 5 pb 1 for e + e~ — > DD* and e + e~ — > DjD* L at 
y/s ~ 10.6 GeV (the subscripts indicate transverse [T] and longitudinal [L] polarization 
of the D*); the cross-section for e + e~ — * DD is expected to be suppressed by a factor of 
~ 10~ 3 . 

This analysis [21] is based on 88.9 fb _1 of data taken at or near the T(45) resonance. 
We reconstruct D° and D + mesons in the decay modes D° — ► K~n + , D° — ► K~n + n + n~ 
and D + — ► K~ii + ti + . D* + mesons are reconstructed in the D°n + decay mode. 




FIGURE 8. Distributions of the mass of the system recoiling against a) D* + , and b) D + . Points with 
error bars show the signal AM reco ii region; hatched histograms correspond to AM reco ii sidebands. The solid 
lines represent the fits described in the text; the dashed lines show the contribution due to events with 
ISR photons of significant energy. The dotted lines show the expected background contribution, c) The 
distribution of M ieco n(D + ) without any requirement on AM reco ii- 



The processes e + e~ — > D^ + D^~ can be identified by energy-momentum balance 
in fully reconstructed events that contain only a pair of charm mesons. However, the 
reconstruction efficiency is small in this case. Taking into account two body kinematics, 
it is sufficient to reconstruct only one of the two charmed mesons in the event to 
identify the processes of interest. We choose the mass of the system recoiling against 
the reconstructed dW (M recoil (Z)M + )) as a discriminating variable: M KCO i\(D^*' + ) = 

^(^fs — E D {*) + ) 2 — P 2 D ^) + , where E D ^)+ and p D (*)+ are the CM energy and momentum 

of the reconstructed . For the signal a peak in the M Teco n distribution around 
the nominal mass of the recoiling D or D* is expected. This method provides a 
significantly higher efficiency, but also a higher background, in comparison to full event 
reconstruction. For the e + e~ — > D^ + D*^ processes we reconstruct in addition a slow 
pion from the D* — > D Q K dow decay. This reduces the background to a negligible level. 

We calculate the difference between the masses of the systems recoil mass against 
a D(*) + 7C~ low combination, and against the DW+ alone, AM recoi i = M recoi i(Z)M + - 
M KC0 [i(D(*) + 7t~ [ow )) . The variable AM reco ii peaks around the nominal D* + — D° mass 
difference with a resolution of OAM recoil ~ 1 MeV as found by Monte Carlo simulation. 
For e + e~ — > D^ + D*~ we combine dW + candidates together with and require 
AMrecoii to be within a ±2 MeV interval around the nominal D* + — D° mass difference. 

The M recoi i(D* + ) and M KCO \\{D + ) distributions are shown in Figs. 8(a) and 8(b), 
respectively. Clear signals are observed in both cases. The higher recoil mass tails in 
the signal distribution are due to initial state radiation (ISR). The hatched histograms 
show the M reco ii distributions for events in the AM reco ii sidebands. 

Since the reconstruction efficiency depends on the production and D*^ 1 helicity an- 
gles, we perform angular analysis before computing cross-sections. A scatter plot of 
the helicity angles for the two Z)*-mesons from e + e~ — > D* + D*~ (cos (D* ec .) vs 
cos0(D* on _ ret .))for the recoil mass region M reco ii(Z)* + ) < 2.1 GeVis shown in Fig. 9(a). 
The distribution is fitted by a sum of three functions corresponding to the D* T D* T , D* T D* L 
and D* L D* L final states, obtained from Monte Carlo simulation. The fit finds the frac- 
tions of DjDj, DjD* L and D* L D* L final states to be (1.5 ±3.6)%, (97.2 ±4.8)% and 
(1.3 ±4.7)%, respectively. Figure 9(b) shows the D* meson helicity distribution for 
e + e — > D + D* . The fraction of the D + D* L final state is found from the fit to be equal 



FIGURE 9. a) The scatter plot cos(0£>* ec ) vs cos(§ D * on _ r J (e + e~ -> D* + D* ). b) £>*+ meson helicity 
angle distribution for (e + e~ — > D + D*~) signal candidates. 

to (95.8 ±5.6)%. 

We search for the process e + e~ — > D + D by studying the recoiling against the 
reconstructed D + (M recoi i). Fig. 8(c) shows the distribution of M recoil (D + ) after D + mass 
sideband subtraction. To extract the e + e~ — > D + D and — > D + D* yields we fit 
this distribution with the sum of two signal functions corresponding to D and D* 
peaks and a threshold function describing background events. The fit finds —13 ±24 
events in the D peak and 935 ±42 in the D* peak. We obtains a e + e~ — > D + D* cross- 
section of 0.61 ± 0.05 pb which agrees with the result using the AM reco ii method. For the 
e + e~ — > D + D cross-section we set an upper limit of 0.04 pb at the 90% confidence 
level. 

In summary, we report the first measurement of the cross-sections for the e + e~ — > 
D*+D* and -> D + D* processes at ^ = 10.6 GeV to be 0.65 ±0.04 ±0.07 pb 
and 0.71 ± 0.05 ± 0.09 pb, respectively, and set an upper limit on the e + e~ — > D + D 
cross-section of 0.04 pb at 90% CL. The measured cross-sections are an order of mag- 
nitude lower than those predicted in Ref. [20], but their relative sizes are as predicted: 
the cross-sections for e + e~ — > D* + D* and e + e~ — > D + D* are found to be close each 
other, while the cross-section for e + e~ — > D + D is much smaller. The helicity decom- 
position for e + e~ — > D* + D* is found to be saturated by the D^D*^ final state (the 
fraction is equal to (97.2 ±4.8)%) and for e+e" -> D+D* — by the D* L final state 
(95.8 ± 5.6%), in good agreement with the predictions of Ref. [20]. 



OBSERVATION OF 7] C (2S) PRODUCTION AND ITS MASS 

MEASUREMENT 

Belle recently observed the t] c (2S) production in exclusive B decays to KKK^K, where 
the T) C (2S) is reconstructed in the K+KgK^ final state. The mass was measured to 
be (3654 ± 6 ± 8) MeV [22] which is much larger than the previous Crystal Ball 
measurement of (3594 ±5) MeV [23]. This year Belle also observed r\ c (2S) production 
(108 ±24 events) in double charmonia events e + e~ — > J/\j/rj c (2S) and confirmed a 



higher T] C (2S) mass [24]. 



CONCLUSION 

We have observed a strong signal for a new charmonium state that decays to k + k~J/ i// 
withM = 3872.0 ±0.6 ±0.5 MeV, Y < 2.3 MeV at 90% CL. We confirm the observation 
of D,y(2317) and D. ? y(2457) and report the first observation of the decay D sJ {2A51) -»• 
D s y. We also observe D s j production in B decays. In B — > tz~ K~ decays all four 
P-wave D** have been observed and their parameters have been measured. For the broad 
DX° and D*,° states there are the first measurements. 
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